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Abstract — The problem of resource allocation is studied for a 
two-user fading orthogonal multiaccess relay channel (MARC) 
where both users (sources) communicate with a destination in 
the presence of a relay. A half-duplex relay is considered that 
transmits on a channel orthogonal to that used by the sources. 
The instantaneous fading state between every transmit-receive 
pair in this network is assumed to be known at both the 
transmitter and receiver. Under an average power constraint 
at each source and the relay, the sum-rate for the achievable 
strategy of decode-and-forward (DF) is maximized over all power 
allocations (policies) at the sources and relay. It is shown that 
the sum-rate maximizing policy exploits the multiuser fading 
diversity to reveal the optimality of opportunistic channel use by 
each user. A geometric interpretation of the optimal power policy 
is also presented. 

I. Introduction 

The multiaccess relay channel (MARC) is a network in 
which several users (source nodes) communicate with a single 
destination in the presence of a relay [1]. The MARC is 
a model for relay-based cooperation in a multiuser network 
where the users have limited power and processing capabilities 
or need tangible incentives to cooperate. We model a MARC 
with a half-duplex relay as an orthogonal MARC where the 
relay transmits on a channel orthogonal to that used by the 
sources (see [2], [3]). The coding strategies developed for the 
relay channel [4] extend readily to the MARC [3]. For exam- 
ple, the strategy of [4, Theorem 1], now often called decode- 
and-forward (DF), has a relay that decodes user messages 
before forwarding them to the destination [5], [6]. Similarly, 
the strategy in [4, Theorem 6], now often called compress- 
and-forward (CF), has the relay quantize its output symbols 
and transmit the resulting quantized bits to the destination [3]. 

We study the problem of resource allocation in a two- 
user ergodic fading orthogonal MARC employing DF under 
the assumption that the instantaneous fading state between 
each transmit-receive pair in this network is known at both 
the transmitter and receiver. Resource allocation for a single- 
user ergodic fading orthogonal relay channel employing DF 
and subject to an average power constraint at the source and 
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relay is studied in [2] (see also [7]). The authors formulate 
the problem as a max-min optimization. They draw parallels 
with the classic minimax optimization in hypothesis testing 
to show that the optimal resource allocation achieves one of 
three solutions depending on the joint fading statistics. The 
orthogonal MARC studied here is a multiaccess generalization 
of the orthogonal relay channel in [2]; however, the optimal 
policies developed in [2] do not extend readily to maximize 
the sum-rate of the MARC. For a two-user MARC, we show 
that the DF sum-rate belongs to one of five disjoint cases 
or lies on the boundary of any two of them. Our results 
reveal two interesting observations: 1) analogously to a classic 
fading multiaccess channel [8], [9], the sum-rate optimal 
policy for each case exploits the multiuser fading diversity 
to opportunistically schedule users; 2) however, these optimal 
policies are not necessarily water-filling solutions. Finally, we 
present a geometric interpretation for each case to highlight 
the effects of node topology in the analysis of multi-terminal 
networks. 

The paper is organized as follows. In Section |II] we model 
the orthogonal MARC with Gaussian noise and fading. In 
Section [III] we present the rate region and determine the power 
policies that maximize the DF sum-rate. Finally, we conclude 
in Section |IV] 

II. Channel Model and Preliminaries 

A two-user MARC consists of two source nodes numbered 
1 and 2, a relay node r , and a destination node d. We write 
K. = {1,2} to denote the set of sources, T = /C U {r} to 
denote the set of transmitters, and V = {r, d} to denote the 
set of receivers. In an orthogonal MARC, the sources transmit 
to the relay and destination on one channel, say channel I, 
while the half-duplex relay transmits to the destination on an 
orthogonal channel 2 as shown in Fig. [T] A fraction 6 of the 
total bandwidth resource is allocated to channel 1 while the 
remaining fraction 6 — 1 — 6* is allocated to channel 2. In the 
fraction 6, the source k transmits the signal Xk while the relay 
and the destination receive Yr and Yd.i respectively. In the 
fraction 9, the relay transmits Xr and the destination receives 
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Fig. 1. A two-user orthogonal MARC. 



Yd.2- In each time symbol (channel use), we then have 

= /i-^ 1^1 ^ hj- 2-^2 ^" (1) 

Yd,i - hd,iXi + hdaX2 + Zd,i (2) 

^d,2 = hd.rXr + Z(i,2 (3) 

where Z^, Zd^i, Zd.2 are independent circularly symmetric 
complex Gaussian noise random variables with zero means 
and unit variances. We write h to denote the vector of fading 
gains, hk.m, for all fc e P and m e T, fc 7^ m, such that h is 
a realization for a given channel use of a jointly stationary and 
ergodic vector fading process iJ. We assume that the fraction 
is fixed a priori and is known at all nodes. 

Over n uses of the channel, the source and relay tranmis- 
sions are constrained in power according to 



for all fc e T. 



(4) 



Since the sources and relay know the fading states of the links 
on which they transmit, they can allocate their transmitted 
signal power according to the channel state information. We 
write Pk (h) to denote the power allocated as a function of the 
channel states h at the k*'^ transmitter, for all k £ T. For an 
ergodic fading channel, (|4]l then simplifies to 



KPkih)) < Pk for all fc e T 



(5) 



where the expectation in (|5]l is over the joint distribution H_. 
We write P (h) to denote a vector of power allocations with 
entries Pk{h) for all k E T, and define V to be the set of 
all P (ft,) whose entries satisfy (|5]l. For ease of notation, we 
henceforth omit the functional dependence of P on h. We use 
the notation C{x) — log(l + x) where the logarithm is to the 
base 2, (x)^ = max(a::, 0), and write Rg = J2kes^k for any 
S QIC. 

III. Sum-Rate Optimal Power Policy 

The DF rate region for a MARC with fixed channel gains 
and a full-duplex relay is developed in [5, Appendix A] (see 
also [6]). For a half-duplex MARC with a fixed h and a 
fixed fraction 9, the DF rate region includes an additional 



conditioning on the half-duplex modes of the relay [3] and 
is the set of all rate pairs (Pi, R2) that satisfy 
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For a stationary and ergodic vector process H_, the channel in 
(HJ-© can be modeled as a set of parallel Gaussian orthogonal 
MARCs, one for each fading instantiation h. For a fixed 
P, the DF rate bounds for this ergodic fading channel are 
obtained by averaging the bounds in ^ and ^ over all 
channel realizations. The DF rate region, TZdf, achieved over 
all P e P, is given by the following theorem. 

Theorem 1: The DF rate region, TZdf, achieved over an 
ergodic fading orthogonal Gaussian MARC is 



Udf = U {7^.(P)n7^d(p)} 

Pev 

where, for all S C IC, we have 



(8) 
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TZr (P) = < (Pi,P2) : Rs<OEC 
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(9) 



and 



nd (P) 



(Pi, P2) : Rs < OEC ^^^^g 



(10) 

Remark 2: The rate region TZd f is convex. This follows 
from the convexity of the set V and the concavity of the log 
function. 

The region TZdf in (Si is a union of the intersections 
of the regions TZr{P_) and TZd{P.) achieved at the relay and 
destination respectively, where the union is over all P_ £ T'. 
Since TZdf is convex, each point on the boundary of TZdf 
is obtained by maximizing the weighted sum /iiPi + IJL2R2 
over all P_ E T', and for all /ii > 0, /i2 > 0. Specifically, 
we determine the optimal policy P* that maximizes the sum- 
rate Pi +P2 when fii — ^2 — ^- Observe from (O that every 
point on the boundary of TZdf results from the intersection of 
T^r{P) and TZd{P) for some P. In Figs. |2]and|3] we illustrate 
the five possible choices for the sum-rate resulting from such 
an intersection. Case 1 and case 2 result when no rate pair on 
the sum-rate plane achieved at one receiver lies within or on 
the boundary of the rate region achieved at the other receiver 
(see Fig. |2]i. On the other hand, cases 3a, 36, and 3c result 
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Fig. 2. Rate region and sum-rate for case 1 and case 2. 
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Fig. 3. Rate region and sum-rate for cases 3a, 3b, and 3c. 



when there is more than one such rate pair as shown in Fig. 
[3j Observe that case 3c corresponds to a boundary case where 
the sum-rate planes overlap. We also consider six boundary 
cases where there is exactly one such rate pair that serves as 
a transition between case 1 or 2 and one of cases 3a, 36 , or 
3c. An example of a boundary case for case 1 and case 3a 
is shown in Fig. |3] We write Bm, ^ V to denote the set of 
P that achieve case i, i = 1,2, 3a, 36, 3c and Bi^n, ^ = 1,2, 
n = 3a, 36, 3c to denote the set of P satisfying each boundary 
case. We show in the sequel that the optimization is simplified 
considerably when the conditions for each case are defined 
such that the sets Bi and Bi^n are disjoint for all i,l,n, and 
thus, are either open or half-open sets such that no two sets 
share a boundary. Finally, we observe that cases 1 and 2 do not 
share a boundary since such a transition (see Fig. |2| requires 
passing through case 3a or 36 or 3c. 

To determine the optimal P*, we first define the conditions 
for each case and determine the policy P^*-' or p(' maxi- 
mizing the sum-rate for case i or the boundary case {l,n). We 
collect the six boundary cases as the last case. The optimal 
policy for each case is determined using Lagrange multipliers 
and the Karush-Kuhn-Tucker (KKT) conditions [10, 5.5.3]. 

Case 1: This case occurs when the power policy P G Si 
achieves the relay and destination regions shown in Fig|2] The 
maximum sum-rate achieved in this case is 



Pi^^ and P^^' as 



(1) 



max (R^f" (P) 
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(11) 



where RJ^"^ 



is the maximum rate achieved by user k at 
receiver m £ I> in (|9]l and (fTOl i. The open set Bi contains 
all P that satisfy 

where P™™ is the rate achieved by user k when it is the 
first user to be successively decoded at a sum-rate corner 
point achieved at receiver m. Since Bi is not known a priori, 
we determine the optimal P*-^-* maximizing Pi + P,2 in ( fTTl i 
over V. Expanding (fTTl l using (|9]l and (fTOl i and applying the 
Lagrange multiplier rule and the KKT conditions, we obtain 



P 



and 



(1) 



{ur^2-wt^y (fc,m)^(l,d),(2,r) (13) 



p(i) 



(14) 



where the water-filling level i^^., fc = 1, 2, r, is determined from 
(|5]l- To ensure that case [T] occurs, we require that P^^^ G Bi, 
i.e., ( flJl l and (fT4l) satisfy (fTSI i. Then, the concavity of the rate 
functions in ^ and ( fTOb suffices to show that P^^^ in ( fT3] l and 
( fT4l ) maximizes ( fTTI ). On the other hand, when P^-^ ^ Si, we 
show that R1+R2 achieves its maximum outside Bi. Note that 
the expression for Pi + R2 for the other cases is not the same 
as that in (fTTT l. The proof follows from the fact that Pi + P2 
in ( fTTT l is a concave function of P for all P €E P. Thus, when 
P'^^ ^ Bi, for every P e Bi there exists a P' e Si with 
a larger sum-rate. Combining this with the fact that the sum- 
rate expressions are continuous while transitioning from one 
case to another at the boundary of the open set Si, ensures 
that the maximum sum-rate is achieved by some P ^ Si. 
We remark that similar arguments also apply to the remaining 
cases, and will be omitted for brevity. Finally, we remark that 
(fT2l i models a network geometry in which the destination and 
source 1 are physically proximal, i.e., they form a cluster, and 
the relay and source 2 form another cluster 

Case 2: The maximum sum-rate achieved for this case is 



where S2 contains all P that satisfy 

RiT (£) < Ri]d (£) and P^"^'^ (P) < P: 
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(15) 



(Z) • (16) 



As in case [T] we can show that ( fTSl ) is maximized by setting 



the optimal P^^"* and Py to the expressions in (fT3] l and (fT4b 
respectively with {k, m) — (1, r), (2, d) provided the resulting 
p(^^ satisfies ( fTSI l. Finally, we remark that the conditions in 
( fT6l l model a network geometry in which the destination and 
source 2 form one cluster while the relay and source 1 form 
another cluster 



3(2) 



Case 3: Consider the cases 3a, 36, and 3c shown in Fig. [3] 
The sum-rate optimization for all three cases simplifies to 



max min ((Rr) 



(17) 



where {Ric)r ™d {Ric)d are the mutual information expres- 
sions in (|9]l and ( fTOl i respectively for S — K., and S3 consists of 
P that do not satisfy ( fT2] l and ( fTSI l either as strict inequalities 
or with equality. We write B3 — B3a U Bsb U Bsc, where Bi, 
i = 3a, 3b, 3c is defined for case i below. The optimization in 
(ITtI i is a multiuser generalization of the single-user max-min 
problem studied in [2]. Extending the known general solution 
to ([TtI i. the optimal policy P^'^^ for this case satisfies one of 
following three conditions 

Case 3a: {Rk)^ IpO) < {RK.)d lp<-*) (18) 
Case 3b: {Ric)r lp(3) > iRic)d IpW) (19) 
Case 3c: {Rfc)r lp(3) = Ipw- (20) 

We proceed to study cases 3a, 3b, and 3c in detail. 

Case 3a: Maximizing {Ric)^ subject to (|5]i results in the KKT 

conditions 

fk<iyM, (21) 



where 



fk — \hr^k\ 




(22) 



For those channel states in which the two users do not see 
the same scaled fading gains at the relay, i.e., |/ir,i|/ t^i 7^ 
\hr.2\/ 1^2, (ISTT l reduces, for all k,j G /C, A; ^ j, and i = 3a, 
to 

, + 





Vk In 2 



(23) 



o.w. 



Observe that the optimal P]^"^^ in 



IS an opportunistic 



water-filling solution that exploits the fading diversity in a 
multiaccess channel with the relay as the receiver. The optimal 
policy at the relay Pr^"-* is given by (O. For P^^a) ^ t^e 
requirement of satisfying (fTsT l. i.e., p(3a) g gg^^ simplifies to 
a threshold condition Pr > Pu (Pi,P2) where Pfe, k £ T, 
is defined in (|5]l and the threshold P^ (Pi, P2) is obtained by 
setting (HI to an equahty. When P^^a) ^ p(3a) ^ ^^^^ 

Pi + P2 is maximized by either case 3b or case 3c. For p(^°) 
^ Z?3, as argued in case [T] the sum-rate is not maximized by 
any P G /B3. Finally, the condition in (fTsl l suggests a geometry 
where the sources and destination are clustered. 
Case 3b : The KKT conditions and optimal P^'^^ for this case 
maximize {RK.)d given by (ISTT i and ( |23] | respectively 

with the subscript in ( |22] | changed to 'd' for aU A: and 
i — 3b. Further, (|2T]) and ( l22b also hold for the relay node, 
k ~ r, and simplifies to the water-filling solution in ( fT4b . 
Thus, as in case 3a, it is optimal to time-duplex the users 
except that this is now based on their scaled fading gains to 
the destination. For P^^*") G S3, satisfying ([191), i-e-. 
G B^b, reduces to satisfying the threshold condition Pr < 
Pi (Pi,P2) where P; (Pi,P2) is determined by setting ( fT9] l 



to an equality. Finally, ( fT9] l implies a geometry in which the 
sources are clustered closer to the relay than to the destination. 
Case 3c (equalizer policy): Maximizing {Ric)^ o'^^r E.^ 
subject to ( |20] | gives P^S'^). This case occurs when p(3c) ^ 
S3c, i.e. P(Pi,P2) < < Pu(Pi,P2)- The resulting 
KKT conditions take the form in (I2TI 1 such that for k ~ 1,2, 



fk 



(1 - a) 



1+ E iv.r^j 



(24) 



For k = r, fr is obtained by replacing in ( |22] | by /id,r 
and scaling by a. The Lagrange multiplier a accounts for the 
boundary condition in ( |20l l and is computed by evaluating ( |20l l 
at p(3c) xjig relay's optimal policy simplifies to the water- 
filling solution in (fT4l i with the first term scaled by a. For 
l^m.il lv\ ^ |/im,2| /i^2, m ~ r 01 d, ( l2Tl i simplifies, for 
in ( l24b and for all fc, j G /C, fc ^ j, and i 3c as 



p(') ^ f (root of /fc|p,=o) fk > fj 



O.W. 



(25) 



From ( |251 ). we see that user k, k — 1,2, transmits opportunis- 
tically over those channel states where a function fk of its 
fading gains and power is larger than that of the other user. 
Note that the optimal policy P^"^ in (|25] | is no longer a water- 
filhng solution at the two sources. Finally, we remark that this 
boundary case occurs for a range of geometries that transition 
from the clustered geometry of case 3a to that of case 3b. 

Case 4: Boundary Cases: We now consider the six bound- 
ary cases and summarize the optimal p('>"), ? = 1, 2, n = 3a, 
36, 3c, for each case. As with case 3c, we observe that the 
geometries for these boundary cases also straddle the clustered 
geometries of the two cases involved. 

Case 1 and Case 3a: The sum-rate optimization for this case 
simplifies to 



max {Ric)r s.t. {Ric) 



omax 

Ri,d 



(26) 



where Si. 3a is the set of all P that satisfy ( fTsT i, (l26T l. and 
one of the inequalities in (fT2] i (see Fig. |3]l. Note that from 
we can write (Pk:)^ R'l™'' + Rf^"^ thus simplifying the 



condition in (|26]l to P""" = P™^'' as shown in Fig. [3] As 
before, we obtain the KKT conditions in ( 1211 1 for k ~ 1,2, 
{k,m) = {l,d), (2,r), and 



fk — 



(1 - a) \hr^k\ 



T'm,k\ 



1+E 



(27) 



3 = i 



where a is the Lagrange multiplier satisfying the boundary 
condition in ( |26] |. When \hrn,i\/i^i 7^ \hrn,2\/i^2, m ~ r or 
d, (I2TJ simplifies p^^'^"-* to the opportunistic non-water-filling 
solution in ( 1251 ) with fk in dZTl i. and for all k, j ^ K., k ^ j, 
{k,m) = (l,d),(2,r), (l,n) = (1,3a). The optimal P^^^'^^^ 
is given by ( fT4] i with the first term scaled by a. 
Case 2 and Case 3a: The sum-rate optimization for this case 
is 



max (P^)^ s.t. P2":;" 



R2T (28) 



where B2.3a is the set of all P that satisfy ( fTSl l, ( |28] l. and 
the remaining inequality in ( fTSl l. The resulting KKT con- 
ditions in (ISTT i use /fc, fc = 1,2, defined in (l27T i but with 
{k,m) = (1, r), (2, d). Note that a captures the equality 
condition in ( |28] ). For 7^ |/im, 21/^^2, m = r or 

rf, and {k,m) — {2,d), p^^'^°^ at each source is given 

by the opportunistic policy in dZSl l. The optimal relay policy 
Pr^''^"'' is the same as that obtained in case 1 and case 3a. 
Case 1 and Case 36: The sum-rate optimization for this case 
is 

max (Pk), s.t. P^T = P^r (29) 

where ;Bi^3b is the set of all P that satisfy (O, (|29]l, and the 
remaining inequality in (fT2] i. The resulting KKT conditions 
satisfy (l2Tl i where /fe, fc = 1,2, is given by (|27] | with 
the subscript 'r' replaced by 'd'. Note that a captures the 
boundary condition in ( l29l l. For |/im,i| ^ |/im,2| /i^2, m — 
r or d, and (k,m) — (2, r), (1, d), p^^'^^-* at each source is 
given by the opportunistic non-water-filling solution in dZSl l. 
The optimal relay policy P^^'^^'' simplifies to the water-filling 
solution in (fl4] |. 

Case 2 and Coie 36: The sum-rate optimization for this case 
simplifies to 

max (Pk), s.t. P™ = P-r (30) 

-rfc02,36 

where 62.36 is the set of all P that satisfy ( fT9b . (l30l l. and 
the remaining inequality in (fTSl l. We remark that the KKT 
conditions are the same as (f2Tl i where fk, k — 1,2, is given 
by (|27]i with 'r' replaced by 'd' and (fc,m) = (l,r),(2,d). 
The resulting p^''^^\ k = 1,2 are given by the opportunistic 
policies in ( fZST l. Finally, the optimal P^^'^'^^ is the same as 
that obtained in case 1 and case 36. 

Case 1 and Ca^e 3c: The sum-rate optimization for this case 
is 

max (Pk), s.t. P^T = R^^"" and P^" = (31) 

-PG/3i,3c ■ ■ ■ ' 

where Si, 3c is the set of all P that satisfy ( f20l ). ( f3T] ). and the 
remaining inequality in (fT2] i. The resulting KKT conditions 
satisfy (f2TI ) for k — 1,2, and (fc. m) = (1, d), (2, r) where 

^ Ifi .|2p. .|2p. 1 I I^TTi.fcl Pk 

1 + E e ' 1 + E e ' + « 

(32) 

and ai, a2, and aa = 1 — ai — a2 are Lagrange multipliers 
that capture the boundary conditions in (fSTT i. When |/im,i| /I'l 
7^ |^m,2| /'^2, m = r or d, for all k, j E K., k j, and pairs 
(fc,m) = (l,d),(2,r), the optimal Pf '^"^^ fc = 1,2, IS given 
by the opportunistic policy in ( fZST i while pj^^'^'^^ is given by 
the water-filling solution in (fT4l i with the first term scaled by 
(ai + 02). 

Cfl.se 2 and Case 3c: The sum-rate optimization for this case 
is 

max {R,c)r s.t. R2T = K'T RT4 = RiT (33) 



where 82,30 is the set of all P that satisfy (|20|, (t33]l, and the 
remaining inequality in (fTST i. The resulting KKT conditions 
and optimal policy p(2.3c) ^j-g j-jjg same as in case 1 and case 
3c but with {k,77i) = (l,r),(2,d). 

Finally, the optimal P* is given by the following theorem. 

Theorem 3: The P* that maximizes the sum-rate is ob- 
tained by computing p(™' or P^^-*-'' starting from case 1 
and proceeding one case at a time, until for some case the 
corresponding p(™) or P^^-'"') satisfies the case conditions. 

IV. Summary and Future Work 

We have developed the power policy that maximizes the 
sum-rate of a two-user orthogonal MARC. We have shown 
that the optimal policy is a function of the channel statistics 
and network geometry and can be classified into two broad 
categories. The first category involves cases where each user 
is clustered with a different receiver as a result of which the 
sum-rate decouples into independent terms for each user and 
the relay. The second category includes the cases where the 
two users are clustered with one of the receivers as well as 
the boundary cases. The first category admits the classic water- 
filling solution at each user and the relay. The optimal policies 
for the second category do not always result in a water-filling 
solution at the sources; however, they reveal the optimality of 
exploiting the multiuser fading diversity to opportunistically 
schedule users. Our results can be generalized to a A'-user 
orthogonal MARC with K > 2. Finally, one could also 
consider the resource allocation problem for CF where the 
challenge lies in solving a non-convex optimization problem. 
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